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1—5 I Abstract 

^ . We compute the 3-point correlation function for a general model of inflation driven by a 

^ ! single, minimally coupled scalar field. Our approach is based on the numerical evaluation of 

1^ \ both the perturbation equations and the integrals which contribute to the 3-point function. 

^ I Consequently, we can analyze models where the potential has a "feature", in the vicinity 

^ ' of which the slow roll parameters may take on large, transient values. This introduces 

^ ■ both scale and shape dependent non-Gaussianities into the primordial perturbations. As an 

'. example of our methodology, we examine the "step" potentials which have been invoked to 

^ I improve the fit to the glitch in the (TT) Ce for £ ~ 30, present in both the one and three 

i-^l year WMAP data sets. We show that for the typical parameter values, the non-Gaussianities 

" " associated with the step are far larger than those in standard slow roll inflation, and may 
even be within reach of a next generation CMB experiment such as Planck. More generally, 

^ ', we use this example to explain that while adding features to potential can improve the flt 

I to the 2-point function, these are generically associated with a greatly enhanced signal at 

I the 3-point level. Moreover, this 3-point signal will have a very nontrivial shape and scale 

^ ' dependence, which is correlated with the form of the 2-point function, and may thus lead to 

■ a consistency check on models of inflation with non-smooth potentials. 
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1 Introduction 



The Cosmic Microwave Background is a treasure trove of information about the primordial 
universe. An exciting set of current and future experiments [H [21 [3] will yield an even richer 
data set, allowing us to make precise tests of inflationary models via their predictions for the 
form of the primordial perturbations. Typically, these predictions are expressed in terms of 
the power spectrum, or 2-point function. For a general set of initial inhomogeneities, further 
information about the perturbations might be gleaned from their higher order correlation 
functions. If the perturbations are exactly Gaussian, then the N-point functions vanish 
exactly when is odd, and is fully specified in terms of the 2-point function for even N. 

All inflationary models predict some level of non-Gaussianity. Moreover, non-linear cor- 
rections will generate non-Gaussianities in the CMB even if the primordial spectrum is 
purely Gaussian Thus any observed non-Gaussianities will be a combination of the 
primordial non-Gaussianity laid down during inflation, and that arising from weak second- 
order couplings between modes. A rough measure of the non-Gaussianities is provided by 
the parameter f^i, 

$(X) = + fNL{^L{xf - ($ (1.1) 

where is the gravitational potential which sources the temperature anisotropics in the 
CMB [5]. Generically, second order couplings between modes yield f^L ~ 0(1), while it 
has been shown in great detail that for standard single field slow roll inflation, /a^l ~ e ~ 
0(10"^), where e is the usual slow roll parameter [B]. The two contributions are cumula- 
tive, and the primordial 3-point function is thus swamped by the evolutionary contribution. 
Moreover, cosmic variance ensures that even a perfect CMB dataset will not allow a conclu- 
sive detection of the 3-point function if f^i is of order unity. In addition, the presence of 
foreground contamination |2l [7j further complicates observational efforts. Current estimates 
suggest that Planck can achieve a limit of I/atlI < 20 ~ 30 |8]. In this paper paper we will 
show that this may permit the detection of interesting features in the 3-point statistics for 
some classes of scale-dependent inflationary models. 

The 3-point correlation function is a separate and independent statistic. Therefore it 
potentially discriminates between models of inflation with degenerate power spectra. More 
optimistically, if the amplitude of the 3-point function is large enough for us to map its 
dependence on both the scale and shape of the momenta triangle, this will be an enor- 
mous boon to early universe cosmology, since the 3-point function encodes information in 
both these properties. However, to realize this possibility, we need to compute the 3-point 
correlation function for a given inflationary model, and compare these predictions to data. 
Here we show how to compute the 3-point function for a general model of inflation driven 
by a single, minimally coupled scalar field. Crucially, we do not assume slow roll, and our 
methodology applies to both simple models with smooth potentials, and more complicated 
scenarios where the potential has an isolated "feature" . Since the 3-point function can have 
both shape and scale dependence, our calculations will underscore the limitations of a single 
scalar statistic such as /nl- Our analysis is purely theoretical, and we do not compare our 
computed 3-point function to data, but this issue is being actively addressed by others [1]. 

As cosmological data improves, we will frequently be faced with the problem of inter- 
preting "glitches" in the observed power spectra, and determining whether these represent 
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genuine departures from some concordance cosmology. For instance, the Ci derived from the 
temperature anisotropics in both the 1 Year and 3 Year WMAP datasets show a noticeable 
departure from the best fit LCDM spectra around I ~ 30, and the fit can be improved by 
adding a small "step" to the inflationary potential. The best-fit parameters for this step were 
extracted from the 1 Year data set in [lOj, and from the 3 Year dataset in |TT], and the two 
analyses find broadly similar values. Of course, any proposal that adds free parameters to 
the inflationary potential is likely to produce a better fit to data. The question then turns to 
whether this improvement is good enough to justify the additional parameters. If we focus 
solely on the 2-point function, this question is addressed via the resulting change in the 
per degree of freedom, or by Bayesian evidence. However, we will see that a potential with 
a sharp, localized feature induces a massive enhancement of the primordial 3-point func- 
tion, relative to our expectations from slow-roll. This enhancement is sufficiently dramatic 
to boost the non-Gaussianities to the point where they may be detected experimentally. 
Consequently, this provides a stringent check on any proposal for a nontrivial inflationary 
potential, since the specific form of the non-Gaussianities will be heavily correlated with the 
2-point function. Such correlations will also extend to higher order statistics such as the 
trispectrum and beyond, though we do not consider them here [111 [131 El ESI IlBl [13 [18] . 

In models with a non-trivial potentials, an accurate computation of the power spec- 
trum must be performed numerically. We begin our examination of the 3-point function 
by reproducing existing semi-analytical results that assume slow roll. We then analyze a 
specific model where the infiaton potential has a small, sharp "step" [19] to illustrate our 
methodology in a non-trivial setting. Since this model has already been used to improve the 
fit between LCDM cosmology and the observed power spectrum, we can take the best-fit 
parameter values for this step and compute the resulting 3-point function. Interestingly, 
although such a feature only causes a small correction in the 2-point correlation function , it 
amplifies the non-Gaussianity by a factor around 1000. Roughly speaking, f^i is boosted 
from (9(.01) to (9(10), making the non-Gaussianity a potentially useful constraint. The 
leading term in the cubic Lagrangian responsible for this boost differs from those used in 
the previous calculations in Ref. [6l [20| [21]. 

The distinctive feature of this non-Gaussianity is its dependence on the scale of the 
momenta triangle [22]. Firstly, it has a characteristic "ringing" behavior; as we will see, 
a suitably defined generalization of /tvl oscillates between positive and negative values. 
Secondly, the impact on the 3-point function is localized around the wavenumbers that are 
most affected by any feature, and dies away over several e-folds. These features distinguish a 
step model from other single field theories with large non-Gaussianities, such as DBI inflation 
[23l [21] or k-inflation models [25], where the non-Gaussianities are present at all scales and 
the running is very slow [26l [22| [2H [271 [28] . 

This paper is divided into the following sections. Section [2| details the step potential 
model that we are considering; we lay out our conventions in this section. Section [3] details 
the analytic forms of the 3-pt correlation functions, which we will proceed to integrate 
numerically in Section [H We discuss our results in Section [5] and conclude in Section [6] 

^As we will discuss later, f^L is not a good statistic for such non-trivial non-Gaussianities, although it 
remains useful as a rough guide to the amplitude of the 3-point function. 
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2 Slow roll model with a feature 



For single field inflation driven by a minimally coupled scalar, the action is 



S 



Ml I ^ 



(2.2) 



where the potential V{(j)) is designed such that the inflaton field (f) is slowly rolling for long 
enough to drive inflation. We define the slow roll parameters bjl^ 



(2.3) 

V = -F^- (2.4) 
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For the quadratic potential, V{(j)) = m'^(f)'^/2, t] = 2e. Assuming transplanckian field values, 
or > Mp, both slow roll parameters and their higher derivative cousins are sub-unity. 
We add a step into the slow roll potential, with the form proposed by 
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(2.5) 



where the step is at 0^, with size c and gradient d respectively. While the presence of the 
step at just the "right" place would require a tuning, one can imagine potentials with many 
steps, so that the odds are high that at least one of them would fall inside the range of (p 
relevant to the cosmological perturbations. 

As discussed in [301 H] this step induces an oscillatory ringing in the power spectrum. 
When the inflaton rolls through the step, it undergoes a strong momentary acceleration. In 
realistic models the step is typically less than 1% of the overall height of the potential. In 
this case, 0^ <^ ^(0) at all times, and e ^ 1. On the other hand, 1] oc V" and its higher 
derivatives can grow dramatically as the inflaton crosses the step, often to the point where 
they exceed unity, as shown in Fig. ([1]). 

We begin with an order of magnitude estimate of the values of the slow roll parameters. 
The step has a depth AV ~ tanh(l)cm202 g^Yid a width A0 ^ 2d (setting Mp = 1). When 
the inflaton falls down the step, AV^ of the potential energy is converted to kinetic energy, 
resulting in an increase in e on the order of Ae ~ AV/H^ ~ 5c, within a time interval 
At ^ A0/0 ^ d/ V cV. Using these quantities, we can estimate 
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^Note these are related to the potential slow roll parameters via ey — (Mp /2){V' /V)"^ , rjv = MpV"/V, 



by e = ey, rj = -2r]i 
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Figure 1: The i]' (left) and e (right) evolution over the step for the model c = 0.0018, d = 
0.022, with its amplitude momentarily t]' ^ lOc^ / [ed'^) . This is the primary source of the 
large non-Gaussianities in the step potential, as the leading term in the 3 point expansion 
is of order ri'e. On the other hand, since the height of the step is small, so the ratio of the 
kinetic energy to the potential energy remains tiny and e <^ 1. 
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(2.7) 



These values have the same order of magnitude as the peaks obtained numerically in Fig. [H 
After this acceleration, the inflaton is damped by the friction term in 



(f> + 3Hcj) + 



dV 
^4 







(2.8) 



and relaxes to the attractor solution. This relaxation time is determined by the first two 
terms in Eq. (12.81) . and is of order H~^. This corresponds to the decay width of the peaks 
in Fig. [U It is easy to see that r] and f] during the relaxation are of order 0{1) and 0{H), 
respectively. Our numerical results were obtained with c = 0.0018, d = 0.022 - the central 
values of the step corresponding to the low-i glitch analyzed by Covi et al. [TT]. Thus 
e ^ 2/02 ~ 2/14. 8^ and get V ^ 10cV(c?^e) ^ 7, consistent with Fig. [H Note that 
aH K, — 1/r is chosen to be around 1 in the plot, r is the conformal time, and prime is the 
derivative with respect to r. 

The non-linear couplings of the perturbations are proportional to powers of the slow roll 
parameters which characterize the background evolution, even when the slow roll parameters 
are large. Consequently, a deviation from slow roll results in a large mixing ("interaction") 
of modes, and large non-Gaussianities. The amplification of the non-Gaussianity relative 
to a smooth potential is roughly characterized by r]feature^t feature/ ^- From our previous 
discussion we know that the rjAt during relaxation after traversing the feature is of order 1. 
During acceleration this quantity will depend on the firm of the feature, and is greater than 
unity for the step potential. Thus any sharp feature will greatly boost the non-Gaussianity. 
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3 3-point correlation functions 



In this section, we sketch out the derivation of the 3-point correlation function, following 
[6l [20l [21], writing in terms of integrals over conformal time r, with the integrand being 
a finite expansion in slow roll parameters. Although we work with the step potential, the 
analysis is applicable to a generic feature. We begin by perturbing the fields in the ADM 
metric i31i 



ds'^ = -N'^dt^ + hij {dx' + N'dt) {dx^ + N^dt) (3.9) 
using the comoving gauge [321 E] 

hij = a^e^'^5ij , (3.10) 

where N and A^* are Lagrangian multipliers, and ( is the scalar perturbation. Note that in 
this gauge 6(f) vanishes. The power spectrum is given by the 2-point correlation function of 
the curvature perturbation 

_ (3.11) 

where is the power spectrum. 

The 3-point correlation function is substantially more complicated. As non-Gaussianities 
arise from departures from nonlinear couplings, we must compute the action (12.21) to cubic 
order in the perturbation: 

S3 = fdtd'x{ah'd' + m^adCf-2mit{dQ{dx) 

+ ^(80(SX)#X + i-mmf + 2/(C)||,} . (3.12) 



where 



X = ahd'''C, (3.13) 
SL . fdd\ 



f{C) = -C^ + terms with derivatives on ( . (3.15) 

Here is the inverse Laplacian and SL/S(\i is the variation of the quadratic action with 
respect to the perturbation (. The cubic action (I3.12p is exact for arbitrary e and rj. Note 
that the highest power of slow-roll parameters which appears is of 0{e^), when we recall 
that X contains a multiplicative factor of e. The fact that the series expansion in slow-roll 
parameters is finite is important since, although e remains small in the presence of a feature. 
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rj and rj' may become large. The last term in (13.121) can be absorbed by a field redefinition 
ofC, 

C-^Cn + /(Cn). (3.16) 

After this redefinition, the interaction Hamiltonian in conformal time is 
HUr) = -ld^x!^ae'a'' + ae'ad(y-2eCmidx) 

+ ^^v'eC + ^(5C)(5X)(5'X) + {^(d'Oidx)'^ . (3.17) 

It is more convenient to work in Fourier space, so 

(C(x)C(x)C(x)> = / ||jg±gL(C(kOC(k,)C(k3)>e--.*-"+-"- (3.18) 

The field redefinition equation fl3.16p introduces some extra terms in the 3-point correlation 
function, 

(C(kOC(k2)C(k3)) = (Cn(kl)Cn(k2)Cn(k3)) 

+ r^(C^(kOCn(k2)Cn(k3)) + sym + 0{v\P^, f) , (3.19) 

where (n{k) denotes the Fourier transform of Cn(^) ^ ^ote that this is not the square of Cn(k). 
In (13.191) . the slow roll parameters are evaluated at the end of inflation. Although we can 
no longer perturbatively expand in terms of ?7 if ?7 is large, the terms that we neglected are 
of higher order in P^. Since ~ 10~^°, we can neglect these unless r] becomes ridiculously 
large, which it never does in any case of interest to us. We have only considered the first 
term in fl3.15p . since all other terms involve at least one derivative of (, and vanish when 
evaluated outside the horizon. 

The 3-point correlation function at some time r after horizon exit is then the vacuum 
expectation value of the three point function in the interaction vacuum 

(C(r,ki)C(r,k2)C(r,k3)) = f dr' a ([C(r, ki)C(r, k2)C(r, kg), ff,„i(r')]) . (3.20) 

Before we proceed, let us pause and point out a subtlety hidden inside equation (13.201) : the 
3-point on the left hand side is evaluated with the interaction vacuum while the the right 
hand side is evaluated at the "true" vacuum. The interaction Hamiltonian Hint evolves 
the true vacuum to the interaction vacuum at the time we evaluate the 3-point function. 
Neglecting this will lead to errors as first pointed out by Maldacena [H [331 El] • 

The three terms in the second line of (13.171) are of higher order in slow-roll parameters, 
and were properly neglected in Refs. [HIED, 21j, which assume that rj and e are always small. 
In our calculation, the terms are small, but the ei]' term may be large and in fact dominant, 
the step. 

We now decompose C, into its Fourier modes and quantize it by writing 

C(r,x) = j (^C(r,k)e^P-, (3.21) 
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with associated operators and mode functions 

C(r, k) = u{t, k)a(k) + u*{t, -k)a^(-k) 



(3.22) 



where a and satisfy the commutation relation [a(k), a^(k')] = (27r)"^5^(k — k'). The "true" 
vacuum is annihilated by the lowering operator a(k). The power spectrum is then 



Pr 



2^ 



(3.23) 



with Mk is M(r, k) evaluated after each mode crosses the horizon. 

Meanwhile v{t, k) is the solution of the linear equation of motion of the quadratic action, 



v'k + k^Vk - —Vk = , 



where we have used the definitions 

Vk = ZUk 



z = a\ple . 



(3.24) 



(3.25) 



Our choice of vacuum implies that the initial condition for the mode function is given by 
the Bunch-Davies vacuum 



Vk{ro) 




2k 



(3.26) 



where we have neglected an irrelevant phase, since (13.241) is rotationally invariant. 

To compute the 3-point correlation function, one simply substitutes these solutions into 
equation fl3.17p . and integrates the mode functions from Tq through to the end of inflation. 
This integral can be done semi-analytically for simple models, provided the slow roll param- 
eters are small and relatively constant. Any departures from this serene picture will render 
the semi-analytic approach intractable - including the step potential we consider here. 

Inspecting equations (13.171) and (13.201) and recalling that Xk oc e(k, it is clear that the 
3-point correlation function consists of a sum of integrals of the form 



L2 oc 3ft 



ler^i OC 3ft 



Wu^{Tend) j rfre?7'a^^i(r)^2(r)^3(7-) 



(3.27) 



(3.28) 



where ^„ is either u^. or du*^ /dr. In a single field model ^^(t) const after Hubble crossing 
as it freezes out, while ^^(t) ^-^kr oscillates rapidly at early times, so its contribution to 
the integral tends to cancel. Thus the integral is dominated by the range of r during which 
the modes leave the horizon. 
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Now let us consider the potential fl2.5l) . We expect e -C 1, making (13.271) negligible. 
However, rj and rj' can clearly be very large as H changes rapidly over a very short time. This 
flagrant violation of slow roll means that we expect a large contribution from the I^rj' term 
in the integral, of the order Arr/'/e x 7^2, where At is the acceleration time for the inflaton 
by the feature. From the qualitative estimation in Sec. [2|, we know this is 7c^/^/((ie^) x 1^2. 
The best parameters from Covi et al. for a step that matches the low-£ glitch seen in current 
CMB data are c = 0.0018 and d = 0.022 Thus we expect a boost in the 3-point statistic 
of (9(100 ~ 1000) at scales affected by the step. Given that the baseline 3-point function is 
proportional to e, we see that this pushes the expected value up to C(10) from (9(10^^). 

We now focus on the J^^/ term, and relegate the discussion of the subleading terms to 
an appendix (where a fleld redeflnition term is also subleading since the effect of the feature 
dies away towards the boundary). So our leading term is 



\ rrend f d \ 

PjMi(T-end) j j draper]' yu*i{r)ul{T)—ul{T) + two permj (27r)^(5^(^ kj) + c.c. 



(3.29) 

where the "two perm" stands for two other terms that are symmetric under permutations 
of the indices 1, 2 and 3, where 1, 2, 3 are shorthand for ki, k2 and k^. 

We must numerically solve the equations of motion to get Uk{T), and then evaluate 
(I3.29p . The contributions from the other flve terms in the interaction Hamiltonian (I3.17P in 
the appendix can be calculated in a similar fashion. 

4 Numerical Method 

The equations of motion for the background flelds and the perturbations are 

H{r) = 3, (4.30) 



da a 



+ 4aV(0)) , (4.31) 



dr^ QM^ 

d'^(f) 2^' (f)' "^^^ (4 32) 

(ir^ a d(j) ' 

where r is our time parameter. We choose initial conditions and units for r such that the 
step occurs around r = — 1, and the mode with co moving wavenumber k = 1 crosses the 
horizon at the same point. The initial conditions of the perturbations are given by (I3.26p . 

We need to evaluate <K7^ and equations f[08l) to f[O0l) . They all have the form iKT!\ 
and fl3.28p . At early times the integrand is highly oscillatory and the net contribution per 
period is thus very small. However the early time limit of the integral must be handled 
with care, since a sharp cutoff imposed by a flnite initial time will introduce a spurious 
contribution of 0(1) into the result. 
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Analytically, the standard procedure is to Wick rotate the integral slightly into the imag- 
inary plane [SI [201 EI] in order to eliminate the oscillatory terms. Unfortunately, it is not 
straightforward to implement this approach numerically. Instead we introduce a "damping" 
factor j3 into the integrand 

1,2 oc f drah^ ^i^2^3 X e^('=i+^'2+'=3)(--re„,)^ (4 34) 

^ '^early 

and similarly for the /^/^ term. This is approximately equivalent to rotating the contour 
integral from the real axis slightly into the imaginary plane r — >■ r(l — ifi). The relative 
error introduced by this numerical procedure can be estimated by taking the difference 
between the integration over e''^h+^^2+k-i)T g(«+/3)(fci+fc2+fc3)T g^^^ jg ^i^us (3. If we ensure 

\(3kTcrossing\ ^ 1 ; (4.35) 
\l3kTearly\ > 1 , (4.36) 

then this contribution is negligible, recalling that r < 0. 

We can adjust both T^ariy and (3. Choosing an appropriate combination will damp out the 
oscillatory contributions at early times while having no effect on the integrand during Hubble 
crossing, which provides the dominant contribution to the 3-point function. We tested this 
scheme against the known analytical results for 3-point functions, which we match to within 
a few percent. However, (3 and Teariy are set by hand in our code, so while this approach is 
suitable for an initial survey, we are developing a more flexible schem^. As an aside, the 
imaginary component of the integral goes to infinity, scaling as a^. However, as long as we 
stop the integration a few efolds after horizon crossing, the integration is stable. 



5 Results and Discussion 
5.1 Shape and Scale 

Every 3-point correlation function has two main attributes: shape and scale. Unlike previous 
treatments, we cannot assume that the 3-point function is scale invariant. Before we compare 
the analytical and numerical results, let us define a more useful parameter than the rather 
unwieldy raw 3-point function, 

(C(ki)C(k2)C(k3)) = (27r)W(ki + k2 + k3)^6;(ki,k2,k3).s (5.37) 

We have factored out the product 1/Hjfcf, since u{t) ~ exp(— zA;r) during inflation and 
there are 6 factors of u in the 3-point function. The function Q depends on the details of 
the integrals listed in Sect. [3] and is some constant that we will come back to later. If 

■^Since finishing this paper, we have developed a more robust numerical methods which regulates the 
integral at early times without using a damping factor |35j . The (3 regularization is prone to suppressing 
the 3-point function since pushing r^ariy to large negative values make the numerical integrations very time 
consuming. In practice the results shown here for the step potential may be ~ 10% lower than their exact 
values, with a greater suppression at larger values of k. 
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the slow roll parameters stay small and constant we can factor the power spectrum out of Q 



(C(ki)C(k2)C(k3)) = (27r)^53(ki + k2 + k3)(Pf)^^^(ki,k2,k3) , (5.38) 

where ^ is a quantity of 0{k^) and the convention here follows Ref. [2T]. In terms of Q, 

NQ = (Pf)U. (5.39) 

Since the power spectrum is statistically isotropic [36], we assume that this is also true 
in the bispectrum and impose rotational symmetry so that, combined with the conserved 
momentum constraint, the 3-point correlation function depends only on the amplitudes of 
the fc's. 

Recall the definition of /atl, equation (II. ip . which is often used when computing the 
non-Gaussianities in the CMB sky |371 |3Hl HI E]- Current constraints on /tvl, from 3 years 
of WMAP data are [39] 

- 36 < /tvl < 100. (5.40) 

As pointed out in [29], this ansatz assumes that the 3-point correlation function has the 
following local form 

(C(ki)C(k2)C(k3)).oca. = (2vr)^53(ki + k, + k3) (^-^f^,{P<f^ ||, (5.41) 
or, in terms of A 

3 

Alocal = - — fNL^ikf. (5.42) 

In other words, the constraint (15.401) is applicable only if the primordial non-Gaussianities 
have the local form fl5.4ip . 

Since equation (11.10 does not contain an explicit scale, this form is scale-invariant (ne- 
glecting the 0{e) correction terms) 

(C(ki)C(k2)C(k3))zocaZ = A^(C(Aki)C(Ak2)C(Ak3)),o,aZ. (5.43) 

Note that there is a small violation of scale-invariance in the 3-point correlation function 
stemming from the power spectrum term which we have ignored in the proceeding argument. 
In fact, equation (15.410 itself hides a further ambiguity: at what value of k are we evaluating 
the power spectrum? This omission is justified if we assume that the power spectrum is 
almost scale-invariant, which is true for most "standard" slow roll single scalar field models. 
However, for the step potential we are considering, the power spectrum itself undergoes 
drastic oscillations of 0{1). We will come back to this point later. 

5.2 Quadratic (m^0^) inflation 

We first reproduce known results for a simple slow roll inflation model [6], 

V{4>) = imV^ (5.44) 
10 



where m = 10 ^Mp. 

Following [B], we can write down the explicit form of A 

AsR = e + ^J^^^jhk'^ + ^^^^^J^^>Jk^k]^ + V (^^^^') + ^(e') ' (5-45) 

where 0{e^) terms include 77'. Comparing this to the local form equation fl5.42p we can see 
that for the equilateral case, — S/atl/IO = lle/8 + 37^/4 = 17e/8 where the second equality 
comes from the fact that rj = 2e for the potential (15.441) . In the other extreme of the 
squeezed triangle case, —S/nl/^^ = 2e + rj = Ae. This result was first derived by Maldacena, 
who pointed out that for single field slow roll inflationary models, /atl = C(e) and is thus 
unobservable. We compared our results to equation fl5.45p . and as Figure [2] shows, our code 
is accurate to within a few percent, which is within the error expected from the analytical 
estimate itself. 

% error 
4 

3.5 

3 

2 . 5 

2 

1.5 

1 

. 5 

1 2 3 4 5 6'' 

Figure 2: Comparison of A/k'^ between the analytical slow roll results of equation fl5.45p 
and numerical results from our code, for the equilateral case where ki = k2 = k^ = k run 
from 0.5 < k < 6.5 and (3 = 0.05. The plot shows the discrepancy between the two sets of 
values. The rapid oscillation at small k indicates the breakdown of the current value of (3 to 
suppress the early time oscillations. Since the k space spans a few efolds, we have included 
the contribution from the e running when computing the analytic estimate but ignored the 
0{e^) terms. As we can see, the two results agree to within a few percent. 

Since the potential has no features, the 3-point correlation function is essentially "scale 
independent" . By this, we mean that triangles which exhibit scaling symmetry posses iden- 
tical values for their 3-point functions. 

5.3 Step potential 

Confident that our code is robust, we now compute the 3-point correlation functions for the 
step potential (12. 5p . We focus on the best fit model of Covi et al. [H] for a step which affects 
modes at large angular scales, namely where (c, 0s, d) = (0.0018, 14.81Mp, 0.022). We choose 
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Figure 3: The power spectrum for the standard m^0^ model and the step potential fl2.5p for 
a decade of k. In this plot, m = 10~^Mp for the m^0^ model while for the step potential we 
have used the parameters (c, 0„ d) = (0.0018, 14.81Mp, 0.022). 



the same inflaton mass m = 10~^Mp used with the m^0^ model, so that we can compare the 
results. 

Consider the power spectrum. Figure ([3]). As first described in [19], the violation of slow 
roll induces a temporary growing/decaying behavior in the evolution of the modes before 
horizon crossing, resulting in a "ringing" of the power spectrum. Thus we cannot simply 
factor out the power spectrum from Q as we have done in equation fl5.38p : doing so will 
introduce spurious contributions to any statistic that we might obtaino On the other hand, 
we expect the 3-point correlation to be of order (-Pf )^, so if we set N = in equation fl5.37p 
where = 4 x 10~^° (i.e. the observed value of the power spectrum), then the dimensionless 
quantity 



53(ki 



P2(27r)' 



■(C(kOC(k2)C(k3)) 



(5.46) 



is what we want to plot. This quantity has the great advantage that the strong k^ running 
from the scaling of the 3-point correlation function (see for example equation fl5.43p ) is 
factored out, so we are left with the scaling that is generated by from the non-linear evolution 
of the modes. 

Consider the equilateral case. Figure (jl]). Analogously to the power spectrum, the pres- 
ence of the step induces a ringing of the equilateral bispectrum. The ringing frequency in 
the 3-point function is roughly 1.5 times that in power spectrum, due to the presence of 
three factors of ( in the integrand instead of two. 



''For example, dividing it by J2i>j PkiPk, [101 SI] will entangle the ringing in the power spectrum with the 
ringing in the bispectrum. The 3-pt correlation function is, in principle, a completely independent statistic 
from the 2-pt correlation function. In our definition, we make it clear that the only content of the division 
by P is simply to state that the 3-pt amplitude is roughly the product of 2-pt amplitude. In the future, we 
envisage that the 3-pt will be part of the set of cosmological parameter space that we can constrain from 
observations of the CMB and large scale structure, and this definition completely disentangle the 2-pt from 
the 3-pt and thus avoid any unphysical cross-correlations between the two statistics. 
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Figure 4: The running of non-Gaussianity Q/kik2k^ (in the equilateral case ki = k2 = 
k^ = k) for the large scale step potential model of {c,(j)s,d) = (0.0018, 14. 81Mp, 0.022). For 
comparison, the standard slow roll model will yield Q/kik2k^ ~ 0{e). 

The plots in Figures illustrate the complicated landscape of the non-Gaussianities. 
Compared to the 0{e) amplitude of the m^0^ model, the non-Gaussianities are magnified 
several hundredfold and can be larger for other viable step parameters (Figure [7]), since we 
are taking the mid-point values from Covi et al. As we explained in Sec. |3l the dominant 
contribution to the non-Gaussianities comes from the t]' term (Figured]). As we estimated, 
Q/kik2k'i is of order 7c^/^/((ie) since the contribution from the 1^2 term is 0{e). This is 
confirmed by our numerical results; the shape of the non-gaussianity remains constant but 
its amplitude is rescaled as we vary c and d together (Figure ([7])). 

In Figure O we plot the integrand (/ci/c2fc3)^-^??'e with respect to r for the modes which 
cross the Hubble horizon around the step. We have suppressed the early time oscillations 
via our "damping trick" in order to isolate the actual contribution to the integral. The 
step also occurs around this time. Comparing this with the rj' plot Figure ([1]), we see 
that the oscillatory nature of u\u2u'^ modulates the single hump of rj' , yielding complicated 
structures. Depending on the combined phases of the modes, the modulation can be both 
constructive and destructive, leading to the "ringing" of the primordial bispectrum. The 
non-linear interactions between the modes have become dominant due to the large coupling 
term erj' around crossing, mixing up the perturbations in non-trivial manner and generating 
large departures from Gaussianity. This is in contrast to the ringing of the power spectrum 
which is caused by the change in the effective mass of the modes as they cross the horizon. 

To bound this specific non-Gaussian signal with current WMAP data, we would have 
to redo the full sky analysis using the data for our complicated scale and shape dependent 
predictions of the bispectrum. This is not a trivial task and requires new methods to be 
developed before we can approach the problem. A brute force computation of the likelihood 
is computationally daunting since the bispectrum here is clearly unfactorizable into a product 
of separate integrals over the A;'^. However, we can make some guesses. We have previously 



recent paper jU has begun to address this important question. 
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Figure 5: The shape of non-Gaussianities Q /kik2k^ for the large scale step potential model 
of {c,(l)s,d) = (0.0018, 14. 81Mp, 0.022), with ki ranging from 6.5 to 2.5 going from top left 
to bottom right. We have set the forbidden triangle regions outside of ki < kj + ki , i ^ j, I 
to —10 for visualization purposes. 
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Figure 6: The real component of the integrand {kik2k';^ya'^eri'ui{0)u2{0)u3{0)ulu2uf^' plotted 
for the modes (/ci, /c2, ^3), with the step occurring around r = — 1 and inflation ending at 
T — 0. Note that the early time oscillations have been suppressed by the inclusion of a 
damping factor P in the integrand. 




-40- -40 



Figure 7: Variation of Q/k^ in response to changing the parameters of the step for the 
equilateral triangle case. The left (right) plot describes the variation of c (d) from the 
smallest to biggest amplitude in the order c = 0.0009, 0.0018, 0.0036, 0.0054 fixing d = 0.022 
{d = 0.044, 0.022, 0.011, 0.0074 fixing c = 0.0018). We keep the location of the step fixed at 
(ps = 14.81Mp. 
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alluded to the fact that the statistic (15.461) is similar to f^i. By comparing equation fl5.42p 
to equation fl5.46p . we see that 

Since the factor 10fcifc2fc3/3Sjfcf is roughly of 0{1) then fp^i ~ Q{ki, k2, ^3)7^1^2^3 ~ C(10) 
which is within sights of the next generation CMB experiments [8j. However, in our case 
only a subset of all possible triangles we could draw on the sky have a large 3-point function, 
which increases the risk that cosmic variance will swamp any signal. On the other hand, we 
have a specific template for the 3-point function, which can then be cross-correlated with 
the 2-point function associated with any putative step. 



6 Conclusion 



While simple models of slow roll inflation do not exhibit observable non-Gaussianities, the 
space of models of inflation that can produce a power spectrum which fits the CMB data is 
infinite. To discriminate between these degenerate models, some other independent statistic 
derived from the CMB must be used and the bispectrum is becoming the most promising 
candidate for this role. Unlike the power spectrum, the computation of the bispectrum for 
a given inflationary model is a messy affair. In this paper, we have developed a robust 
numerical code which can compute the 3-point correlation function for a non-slow roll single 
field model. 

Using this method, we numerically integrated the 3-point correlation functions for a 
class of inflationary models where slow roll is violated for a brief moment. We show that 
the addition of this step breaks the scale invariance of the bispectrum in a non-trivial way, 
leading to large non-Gaussianities. This numerical technique can be extended to solve more 
complicated problems such as multifield models of inflation [121 112] and to compute higher 
statistics such as the trispectrum [16]. 

Mirroring the behavior of the power spectrum, a temporary violation of slow roll in- 
troduces a ringing into the bispectrum, resulting in a structure much like the transient 
vibrations of a rug being beaten. These large non-Gaussianities are generated when the 
non-linear coupling between the modes, erj', becomes temporarily large during this period of 
non-slow roll behavior. In order to check whether we can constrain such models with data, 
a full sky analysis must be undertaken with either current or future data, a task which is 
beyond the scope of this paper. However, we argue that at least for the models considered 
in this paper, the non-Gaussianities may be large enough to be within reach of the Planck 
satellite. 

In the near future we can expect even better quality CMB data to become available, 
opening up the possibility of cross-correlating higher order statistics such as the bispectrum 
with the power spectrum. This should provide a powerful consistency check on any models 
of inflation with a sharp feature in the potential, even if this improves the fit to the 2-point 
function. Consequently, we have laid the groundwork for an analog of the "consistency 
relationship" of simple models of slow roll inflation, where the amplitude and slope of the 
scalar and tensor spectra are expressed as functions of just three independent parameters - 
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the observed amplitude of the scalar spectrum, and the first two slow roll parameters, e and rj. 
In this case considered here, the parameters that describe any step or other near-discontinuity 
would fix both the 2- and 3-point functions, yielding a consistency test for a very broad class 
of inflationary models. In practice we could fit separately to the 2- and 3-point functions 
and then check that the results matched. Alternatively, one could simultaneously include 
information from both the power spectrum and bispectrum into the cosmological parameter 
estimation process. In this case, we could fit directly to the parameters c, d and (ps (along 
with the "average" slow roll parameters) , generalizing the approach of ^51 H6] . Finally, if 
it ever becomes possible to determine the 3-point function from observations of large scale 
structure, this information would further tighten the constraints on models where slow roll 
is briefly violated. 
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A Other terms 

In Section [3], we have given the expressions for the leading non-Gaussianity in the case when 
7] is large but e remains small. In this appendix, we give the expressions for the subleading 
terms for this case from the Hamiltonian (13.171) . 

The subleading contributions come from the first line of (13.171) . Contribution from a^e^CC^ 
term: 

/fend / -I— r \ / du* du* \ 

dr ah^ \ T[ Ui{Tend)j \ + two permj {2Tif5\Y^ k,) + c.c. . (A.48) 

Contribution from ae^C,{dCf term: 

/Tend [ \ 

dr ah^ \ \{ui{T,nd)u*{T)\ (ki • ^2 + two perm) {2iTf6\^\ii) + c.c. . (A.49) 

Contribution from —2atC,{dC){dx) term: 

(A.50) 

The redefinition fl3.15p contributes 



-2i 



|2l, |2 



2 l^2| IM3I 



{2Tif5^(Y^ ki) + two perm . (A.51) 
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Note that in Ref. [SI [201 EI], where sharp features are absent, these four terms give leading 
contributions to non-Gaussianities. 

The last two terms in fl3.17p are further suppressed by e, but for completeness we list 
their contribution in the following. Contribution from [e/2a)d(dxd'^x term: 



(A.52) 



2 

Contribution from (e/4a)(9^C)((9x)^ term: 



(A.53) 
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